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Abstract.  Scattering  of  Rayleigh-Lamb  waves  by  a  normal  surface-breaking  crack  in  a 
plate  has  been  studied  both  theoretically  and  experimentally.  The  two-dimensionality  of 
the  far  field,  generated  by  a  ball  impact  source,  is  exploited  to  characterize  the  source 
function  using  a  direct  integration  technique.  The  scattering  of  waves  generated  by  this 
impact  source  by  the  crack  is  subsequently  solved  by  employing  a  Green’s  function 
integral  expression  for  the  scattered  field  coupled  with  a  finite  element  representation  of 
the  near  field.  It  is  shown  that  theoretical  results  of  plate  response,  both  in  frequency 
and  time,  are  similar  to  those  obtained  experimentally.  Additionally,  implications  for 
practical  applications  are  discussed.  v  f'S  .  >  y 

f  distribution  stv^ : - 
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Introduction  Dumpur.;-  -  ■> 

Propagation  of  guided  Rayleigh-Lamb  waves  in  a  plate  is  of  interest  in  seismol¬ 
ogy,  electrical  devices,  ultrasonic  material  characterization,  and  ultrasonic 
nondestructive  evaluation  of  defects.  There  have  been  numerous  investigations 
of  this  problem  since  the  early  works  of  Rayleigh  and  Lamb,  and  reviews  of  the 
early  literature  can  be  found  [1,  2],  In  recent  years  attention  has  been  focused 
on  the  surface  response  of  a  plate  due  to  buried  or  surface  sources  in  the 
context  of  acoustic  emission.  The  direct  problem  of  the  response  of  a  plate  due 
to  various  sources  has  been  studied  by  many  authors  [3-10]  for  the  purpose  of 
analyzing  signals  from  acoustic  emissions.  Experimental  measurements  of  sur¬ 
face  response  due  to  normal  force  sources  acting  on  a  plate  have  also  been 
examined  [10-16],  In  this  paper  we  have  theoretically  studied  the  surface  re¬ 
sponse  of  a  plate  due  to  a  vertical  transient  line  force  acting  on  the  surface  of 
the  plate.  The  problem  is  considered  to  be  two-dimensional,  i.e.,  plane  strain. 
The  model  studied  is  a  simulation  of  the  £*pvtf>ental  set-up  shown  in  Fig.  1 . 
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Fig.  1.  The  glass  specimen  with  support  structure  and  ball  drop  source  apparatus.  The  structure 
supports  the  25.4  mm  by  5.6  mm  by  2.28  m  specimen  but  is  mechanically  decoupled  by  a  high 
attenuation,  impedance  mismatch  urethane  foam.  This  minimizes  the  loss  of  energy  to  the  support¬ 
ing  structure  while  isolating  the  specimen  from  external  vibrations.  Aluminum  coated  microscope 
slides  are  glued  to  the  sample  surface  to  provide  the  electrically  conductive  and  smooth  surface 
required  by  the  sensor. 


Theoretical  Formulation  and  Solution 

The  response  of  an  elastic  plate  due  to  an  appied  source  of  excitation  is  charac¬ 
terized  by  the  type  of  source.  Theoretical  knowledge  of  response  due  to  known 
sources  helps  in  identifying  source  functions  generated  in  experimental  set-ups, 
or  in  practical  applications.  Characterizing  a  physical  source  in  a  numerical 
form  is  a  necessary  step  towards  analyzing  the  scattered  field  for  nondestruc¬ 
tive  evaluation  of  defects  in  a  plate. 

The  propagation  of  waves  in  elastic  plates  has  long  been  studied  by  many 
authors.  Mindlin  [1]  gave  a  thorough  discussion  of  the  solutions  to  the  govern¬ 
ing  equations  of  wave  propagation  in  isotropic  plates.  Since  then,  many  studies 
have  been  reported  on  the  forced  motion  of  a  plate.  Recently  transient  wave 
propagation  in  an  infinite  plate  due  to  a  point  source  was  studied  by  Ceranoglu 
and  Pao  [2].  They  employed  a  generalized  ray  theory  to  compute  the  desired 
response;  however,  this  method  is  very  cumbersome  when  considering  re¬ 
sponses  in  the  far  field,  i.e.,  several  plate  thicknesses  away  from  the  source.  To 
avoid  this  difficulty  Weaver  and  Pao  [8]  used  a  normal  mode  analysis  to  analyze 
responses  in  the  far  field.  The  response  of  a  finite  plate  was  described  by  a 
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double  summation  over  wave  number  and  eigenfrequencies.  Letting  the  radius 
approach  infinity,  the  summation  over  wave  number  gets  replaced  by  its  indefi¬ 
nite  integral.  Vasudevan  and  Mai  [9]  presented  a  method  based  on  the  classical 
integral  transform  technique.  The  response  at  a  particular  frequency  is  ob¬ 
tained  by  using  the  method  of  residues  at  the  poles  (real,  imaginary,  and  com¬ 
plex)  which  represent  the  modes  of  propagation  discussed  in  [I].  Finally,  an 
inverse  fast-Fourier-transform  (FFT)  is  used  to  obtain  the  time  response  of  the 
plate. 

In  this  work  we  consider  the  two-dimensional  problem  of  a  uniform  line 
force  acting  on  one  fact  of  a  plate.  The  response  of  the  plate  due  to  an  impulsive 
load  is  written  in  the  form  of  a  double  integral  over  the  frequency  and  wave 
number.  For  any  frequency,  direct  integration  is  carried  out  with  respect  to  the 
real  wave  number  after  adding  a  small  imaginary  part  to  the  real  frequency.  The 
Fourier  spectrum  of  the  experimentally  generated  source  function  is  then  ob¬ 
tained  by  calculating  the  ratio  of  the  observed  response  and  the  input  response. 
The  corresponding  representation  in  time  domain  is  then  carried  out  by  taking 
the  inverse  FFT  of  the  frequency  response.  This  numerical  representation  of 
the  source  function  is  used  to  study  the  scattering  of  waves  by  a  normal  surface 
crack  in  a  plate. 


Evaluation  of  Green’s  Function 

The  Green’s  function  for  a  medium,  G,y(x,  z,  f;  x',  z',  0),  is  defined  as  the 
displacement  in  the  iIh  direction  produced  at  a  point  (jc,  z)  at  time  t  by  an 
impulsive  line  force  of  unit  magnitude  acting  at  t  =  0  in  the  jth  direction.  For 
computing  Gy  we  consider  the  source  to  be  uniform  along  the  y-axis,  i.e. ,  we 
assume  the  problem  to  be  a  2-D  plane  strain  problem.  In  index  notation  the 
Green’s  displacement  vector  Gj  must  satisfy  the  equation  of  motion: 

mV2G j  +  (X  +  p)\V  ■  G;  =  pGj  -  e,8(f)8(jt  -  x')8(z  -  z'),j  =1,3.  (1) 

Assuming  the  displacement  vector  to  be  of  the  form 

Gy  =  Vd>y  +  VA(^e2),  (2) 


one  gets 


_  d<f)}  di hj 

G]J  -  *7  ~  IF 


_  d<t>j  dibj 

Gij  "  aF  a*’ 


where  <J>;  and  t|/y  satisfy  the  differential  equations  (except  at  x\  z ' ): 

V2<bj  =  $/  •  =  ^2  $/' 

C?  =  (A  +  2 ix)/p.  Cl  =  pip. 


(3) 
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Fig.  2.  Two-dimensional 
layered  plate  with  top  and 
bottom  surfaces  traction  free. 
The  layered  plate  contains  n 
plates  of  thickness  h ,  welded 
together  with  total  thickness 
of  nh  =  H. 


Free  Surface  x 


For  generality,  let  us  consider  n  plates  of  equal  thickness  h  that  are  welded 
together  with  the  top  and  bottom  surfaces  being  traction  free  (See  Fig.  2).  In 
order  to  solve  for  the  Green’s  displacement,  Gy,  in  the  n-layered  plate  we  take 
the  time  dependence  first  in  the  form  where  tv  is  circular  frequency. 

Furthermore,  we  write 


<t>y  =  J  d>y  eikx  dx 

%  =  J\  %eikx  dx.  (4) 

Then  it  is  easily  shown  that  within  the  mlh  layer  d>y  and  'Fy  can  be  written  as 

<t>y"  =  ATeir""z  +  Bj’e1 +  5™  (5) 

'F™  =  Aj'e‘^z  +  BJ'  eirnJh~z)  +  S",  (6) 

where 

V\m  =  -  k2,  r)2m  =  Vk22m  -  k2, 

k\m  wlC\ m ,  k2m  wlC2 m 

S']'  and  T”  depend  on  the  direction  and  location  of  the  applied  unit  load.  For  a 
force  in  the  z-direction  acting  at  (x1 ,  z')  they  are  described  as 


5^  =  ±^imehim\z-z'\ 

j-™  —  Bime‘Vim^z~z^ 


(7) 


Scattering  of  Rayleigh-Lamb  Waves 


201 


where 


p-UiX'  K  p-W 

Aim  —  _  Dim  _ _ 

(2pmfc|J’ 

In  the  above  the  ±  signs  correspond  to  z  >  or  <  z\  respectively.  Similarly,  for 
a  force  in  the  x-direction  at  (x\  z')  we  get 


Tf  =  ±Blmeir,2jn  l£_z’l  (8) 


with 


Aim  =  “  _ ; _  Dim  —  .  . 

TJlw  (2flmk\my  (2 IAmk\my 

If  there  is  no  applied  load  in  the  mth  layer,  then  Sf  =  Tf  =  0.  Note  that  in 
writing  the  above  solution  a  local  coordinate  system  has  been  chosen  with  the 
origin  on  the  top  of  the  layer. 

The  stresses  associated  with  the  Green's  displacement  G *  are  calculated 
from  the  stress-strain  relations.  The  relevant  stress  components  needed  to 
satisfy  the  continuity  conditions  at  the  interface  and  the  free-surface  conditions 
are 


XU  -  (I  c;  f  l  c3) 

SSf  =  ■  G7  +  2fim  j:  G”.  (9) 

d<. 

Using  Eqs.  (5)  and  (6)  in  Eq.  (3),  and  in  Eq.  (9),  the  Fourier  transforms  of 
the  displacement-stress  vector  in  the  mth  layer  can  be  written  in  matrix  form  as 

[Ujm(z)1  =  [lm(z)][X?]  +  [  Y?\,  (10) 

where 


[U?(z) ]  =  [G7j,  G",  ±Zt,  2”]r 

[Xj]  =  [A?C?B?D?]T 

Expressions  for  lm  and  Yf  are  given  in  Appendix  I.  The / represents  the  Fourier 
transform  of  fix). 

The  boundary  conditions  to  be  satisfied  by  [£/;]  are 
%lZZj  =  =  0  at  z  =  0 


(ID 
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[U?(h))  =  [t/j"+l(0)],  m  =  l,  ....  n  -  1  (12) 

=  ±b  =  0  at  z  =  h  (13) 

These  give  4 n  conditions  for  the  4 n  unknown  constants  [A 'f],  m  =  1,  n. 

Assembling  the  Eqs.  (11 )— ( 13)  into  a  global  matrix  we  obtain  the  matrix 
equation 

l/-’][A]  =  [T]  (14) 


where  [F]  is  a  4 «  x  4 n  matrix, 

[X]  =  ...,  A"  ,C"  ,B"  ,Dj]T  and  [  T]  =  [IT, .  YYn+l)T 


with 

YY,  =  [  —  Ky  (3),  —  Tj(4)]j=0 
YYn+\  =  [- 17(3),-  T;(4)][=A 
YYm  =  [17Uo- Tr'U]r,  m  =  2 . «. 

Note  that  YY,  and  KT„+i  are  2  x  1  matrices  and  the  rest  are  4  x  1. 

Solving  Eq.  (14)  for  [A]  we  obtain  ail  of  the  constants.  Once  these  con¬ 
stants  are  known  the  displacements  and  stresses  at  any  point  can  be  found.  It 
may  be  noted  that  setting  the  determinant  of  [F]  to  zero  gives  the  dispersion 
equation  that  relates  the  frequency  k2\  to  the  wave  number  k.  This  has  been 
discussed  recently  by  Mai  [17]. 


Extraction  of  the  Source  Function 

We  shall  now  focus  our  attention  on  a  single  homogeneous  plate  and  obtain  an 
expression  for  the  z-component  of  the  displacement  at  z  =  0  for  a  vertical  load 
also  acting  at  z  =  0.  This  will  be  denoted  by  D(x ,  x\  f)  for  a  particular 
frequency  w  =  2irf,  /being  measured  in  cycles  per  second.  Following  the 
procedure  for  an  n-layered  plate  we  get  the  four  constants  A\ ,  Cj,  and  D\. 
Using  these  constants  we  find  from  Eqs.  (5),  (6)  and  (3)  (see  also  Achenbach 
(18]): 

D{X'  X  ';f)  =  2^  /I  {l2'A'3  +  l'22C'3  +  12383  +  l23°3  +  rliJU-o**  dk  (15) 


where  A3,  ...,  D\  are  the  solutions  to  the  equation 


a 

C  ae‘Vu>< 

'M' 

aA3]  +  cB 31 

b 

—  Q  —fjgiVllh 

3 

c\ 

►  =  i 

~bA3'  +  aB3' 

ceIT|JI  *  a 

—c 

aA3  *eir,"h  4-  cB3le'Vl,h 

beiy)"h 

—ae" if*  -b 

~~a 

Pi 

-bA3xeiru'h  +  aB3>e'm'hJ 
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The  integration  is  done  numerically  according  to  an  adaptive  scheme  that  is 
a  modification  of  one  by  Xu  and  Mai  [19].  Using  this  scheme  we  evaluate 
displacements  and  stresses  as  a  function  of  k2  at  any  point  (jc,  y)  for  a  source  at 
(x’,  z')  in  the  plate.  The  results  of  integration  were  checked  for  some  arbitrary 
frequencies  with  those  obtained  by  the  method  of  modal  summation.  In  this 
particular  problem,  since  the  source  and  receiver  are  at  the  top  surface,  z  and  z' 
are  zero. 

Experimentally  we  measure  the  time  responses  of  normal  surface  displace¬ 
ments  at  the  receivers  r\  and  r2  that  are  distances  d\  and  d2  from  the  source, 
respectively.  The  source  is  considered  to  be  a  line  source  S  uniform  along  the 
width  of  the  plate  and  acting  in  the  z-direction.  We  then  obtain  the  correspond¬ 
ing  frequency  responses  /?*(/)  and  R2(f)  by  means  of  an  FFT.  Now,  the  re¬ 
sponse  at  a  frequency/,  R(f)  at  the  point  (x,  0)  is  given  by  S(/)*D(/),  where 
£(/)  is  the  frequency  response  of  the  source  and  />(/)  is  given  by  Eq.  (15).  We 
obtain  D(f)  analytically  and  hence,  compute  S(f)  from  the  observed  /?(/). 

Thus  we  get  S  for  discrete  frequencies /=  n  x  A/,  where  n  =  1,2 . m,  m  = 

2’,  I  an  integer,  and  A /  =  1/7},  7}  being  the  length  of  the  time  window.  This 
S{f),  which  is  a  complex  number,  represents  the  source  function  in  the  fre¬ 
quency  domain.  The  corresponding  time  representation  of  the  source  is  ob¬ 
tained  by  taking  the  inverse  FFT.  For  two  different  observation  points  the  time 
signals  measured  by  the  receivers  at  ( x ,,  0)  and  (x2,  0)  are  labeled  as  7,  and  72, 
respectively,  with  their  corresponding  frequency  spectra  as  /?,(/)  and  R2(f). 
From  the  theoretically  computed  D\(f)  and  D2(f)  at  these  points,  since 

W)  =  Rt(f)IDt(f),  S2(f)  -  R-JJ)ID2(f), 

we  thereby  obtain  the  source  function  from  observations  at  different  receiver 
positions.  Since  the  source  remains  the  same  for  both  the  receiver  positions, 
S|(/)  and  S2(f)  should  be  identical,  at  least  in  theory.  In  practice,  they  are 
found  to  agree  extremely  well  when  the  receivers  are  more  than  a  few  plate 
thicknesses  away,  (i.e.,  more  than  3H).  Note  that  we  did  not  evaluate  the 
source  functions  at  zero  frequency,  so  the  corresponding  time  representation  of 
the  source  will  have  an  inherent  DC  component. 

The  positions  of  the  receivers  considered  in  this  work  are  x,  =  5H  and  x2  = 
8H.  Since  the  time  signals  do  not  decay  sufficiently  at  the  end  of  the  time 
window,  a  window  is  used  on  the  experimental  time  signals  to  force  them  to 

decay.  The  window  used  is  cos(/  n/(2 n)),  where  i  =  1.2 . n,  n  was  taken  to 

be  512.  This  window  is  imposed  on  all  of  the  experimental  time  signals  for 
computational  purposes.  Additionally,  zero’s  were  added  to  the  time  signals  to 
increase  the  length  of  the  signals  by  a  factor  of  2.  This  introduction  of  zero's 
enables  the  A/to  be  reduced  by  a  factor  of  2,  for  increased  numerical  accuracy, 
without  sacrificing  the  integrity  of  the  information  contained  in  the  signal. 


Scattering  by  a  Surface-Breaking  Crack 

The  problem  of  scattering  of  incident  body  waves  by  a  surface-breaking  crack 
is  solved  by  employing  a  hybrid  method  which  combines  the  advantages  of  the 
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Fig.  3.  Regions  and  boundaries  used 
by  the  hybrid  method  which  com¬ 
bines  the  advantages  of  the  finite 
element  technique  and  the  boundary 
integral  method.  The  boundaries  are 
B  and  C,  with  C'  being  the  free 
surface.  The  regions  are  R ,  (inte¬ 
rior)  and  R0  (exterior).  The  finite 
element  mesh  that  is  contained  in  /?, 
is  illustrated  in  Fig.  7. 


finite  element  technique  and  the  boundary  integral  method.  The  details  of  the 
hybrid  method  are  discussed  by  Khair  et  al.  [20,  21],  who  considered  scattering 
in  a  semi-infinite  medium.  For  this  purpose  we  consider  two  artificial  bound¬ 
aries  C  and  B  (Fig.  3).  The  medium  is  now  divided  into  two  regions.  The 
interior  region  R ,  is  bounded  by  B,  and  the  free  surface.  The  exterior  region  R0 
is  bounded  by  the  free  surface,  the  boundary  C,  the  boundary  z  =  H  and 
extends  to  infinity  in  the  ±x  directions.  The  area  between  C  and  B  is  shared  by 
both  regions.  In  regions  R0  and  R i  the  governing  equation  of  elastic  motion  is 
written  as,  assuming  the  time  dependence 

Tijj  +  pw2Ui  =  0,  /,  f  =  1,3  (16) 

where  Ty  is  the  stress  tensor,  p  the  mass  density  and  (/,  the  i,h  displacement 
component.  Solution  to  Eq.  (16)  satisfying  the  stress-free  boundary  conditions 
along  the  surface  of  the  half-space  and  the  crack  surface  is  sought. 


Solution  in  the  Exterior  Region  R0 

In  this  region  the  displacement  is  composed  of  two  parts. 

Ui  =  Uf  +  U\s)  (17) 

where  £/i0,(/  =  1,3)  represents  the  free  field  displacement  components  (the 
incident  field  due  to  the  line  load)  and  V)s)  is  the  scattered  field.  The  scattered 
displacement  field  in  R0  is  represented  by  a  surface  integral  [20],  after  dropping 
the  factor  e~iwl. 


Uy\x\  z')  =  fc  (Gan  -  U,2,jk)nk  dS  (18) 

Here  Gy  is  the  previously  described  Green’s  displacement  tensor  and  nk  defines 
the  components  of  the  outward  unit  normal  vector  to  C.  The  integration  along 
C  is  carried  out  in  the  clockwise  direction. 
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Solution  in  the  Interior  Region  Ri 

This  region  is  divided  into  finite  elements  having  /V/  number  of  interior  nodes 
and  Nb  number  of  boundary  nodes.  For  the  finite  element  representation  in 
region  /?,  the  energy  functional  is  taken  to  be 

E  =  1/2  f  f  [T  •  €*  -  pw2U  •  U*]dA  -  1/2  [  [tfl  •  US  +  t B-  UB]  dS  (19) 

J Jr  i  jb 

where  *  denotes  complex  conjugate  and  T,  e  are  stress  and  strain  vectors 
defined  as 


T=  [Txx,  Tzz,  Txl]T  (20) 

e  =  [e«,  ez;,  e^]7.  (21) 

Superscript  T  denotes  transpose.  tB  and  UB  denote  the  traction  and  displace¬ 
ment  at  contour  B,  respectively.  It  is  assumed  that  the  displacement  field  within 
an  element  is  represented  in  terms  of  the  shape  functions  z)  and  elemental 
nodal  displacements  Utf] as 


£/[f) 


1,  3) 


(22) 


The  number  of  nodes  in  each  element  is  given  by  n.  Substituting  Eq.  (22)  into 
the  strain-displacement  relations  and  these,  in  turn  into  the  stress-strain  rela¬ 
tions  provides  T\f  and  ejf1.  Substituting  these  in  Eq.  (19)  and  taking  its  varia¬ 
tion,  the  equation  of  motion  for  region  R i  can  be  written  as 


Sii 

Sib 

fu-l 

l-{0] 

_Sbi 

Sbb_ 

luj 

IYbJ 

The  elemental  impedance  matrix  [5]P  is  given  by 

[Sir  =  (J  {[ Bt)T[D][Be ]  -  pw2[<t>e]T{<t>e]}  dA  (24) 

and  Yg  is  the  nodal  force  vector  due  to  surface  tractions  on  the  boundary. 
Using  the  top  set  of  equations  in  Eq.  (23)  we  get 

Wt]  =  -[S„]  '[S/bKUb}  (25) 

Matrices  [£,]  and  [D ]  have  been  derived  before  [20].  Combining  Eqs.  (18)  and 
(25)  the  boundary  displacement  UB  is  found  to  be 

{UB}  =  (-[Ab/][S„]-'[S/b]  +  AbbMIM  +  {Uf}.  (26) 

Once  {Ub}  is  found  by  solving  Eq.  (26),  { U /}  is  found  from  Eq.  (25).  The 
aforementioned  method  was  used  to  solve  for  the  Z-component  of  the  surface 
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displacements  in  the  presence  of  a  normal  surface-breaking  crack  when  the 
plate  is  excited  by  a  source  on  the  same  surface  as  the  crack.  The  parameters 
used  for  Poisson’s  ratio.  Young’s  modulus,  and  the  P  and  SV  velocities  were 
those  of  standard  window  pane  glass  (see  the  following  Experimental  Proce¬ 
dure). 


Experimental  Procedure 

In  the  experimental  part  of  this  study  we  generate  Rayleigh-Lamb  wa''es  in 
glass  samples.  Our  purpose  is  two-fold.  On  the  one  hand  we  wish  to  create  an 
approximation  to  the  two-dimensional  case  which  is  theoretically  tractable. 
Here  we  need  to  determine  the  source  function  and  study  the  effect  of  a  crack 
on  the  propagation  of  the  waves  in  the  far  field.  On  the  other  hand  we  are 
exploring  other  possible  experimental  arrangements  with  the  aim  of  producing 
a  practical  tool  for  crack  detection.  In  these  latter  arrangements  we  do  not  limit 
ourselves  to  those  configurations  which  lend  themselves  to  theoretical  treat¬ 
ments. 

We  will  now  describe  the  experimental  apparatus  which  consists  of  the 
following  parts:  the  source,  the  receiver,  the  sample,  and  the  data  acquisition 
system. 


The  Source 

The  source  consists  of  1.6  mm  steel  ball  bearings  impinging  on  the  glass  sur¬ 
face.  The  balls  are  dropped  through  a  carefully  reamed  hole  in  a  brass  tube.  The 
hole  diameter  is  just  slightly  larger  than  the  ball  diameter  to  assure  laminar  air 
flow  around  the  balls  and  thus  reproducible  ball  impact  velocities.  The  balls  are 
held  by  a  magnet  and  released  by  removing  the  magnetic  field  (Fig.  1). 


The  Receiver 

Since  we  are  comparing  normal  surface  displacements  obtained  from  experi¬ 
ments  with  theoretically  derived  results,  a  transducer  (receiver)  with  a  known 
transfer  function  is  required,  i.e. ,  one  with  an  absolute  calibration  of  amplitude 
response  versus  frequency.  We  have  developed  such  a  transducer  (see  [22]), 
that  alleviates  many  of  the  problems  which  exist  with  conventional  piezoelec¬ 
tric  transducers.  Electronically  this  capacitive  transducer  (CT)  is  based  on  the 
VideoDisk  player  development  by  RCA  (Radio  Corporation  of  America).  It 
detects  the  changes  in  a  100  nm  gap  between  a  0.3  mm  diameter  needle  and  the 
surface  of  the  sample.  It  responds  only  to  the  component  of  the  signal  that  is 
normal  to  the  surface.  The  frequency  response  is  flat  within  7  dB  from  below  10 
kHz  to  above  6  MHz.  The  sensitivity  of  the  instrument  approaches  that  of 
piezoelectric  transducers  at  0.44  V/nm  with  a  noise  level  of  0.004  volts.  At  a 
one  to  one  signal  to  noise  level,  its  absolute  sensitivity  is  therefore  about  0.015 
nm.  The  dynamic  range  covers  over  three  orders  of  magnitude  from  0.015  nm 
to  over  50  nm.  The  sensitivity  of  the  CT  is  given  in  Fig.  4  as  a  function  of 
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Frequency  [kHz] 

Fig.  4.  Sensitivity  of  the  CJ  versus  frequency.  Surface  displacements  were  generated  by  a  resonat¬ 
ing  rod  excited  by  a  piezoelectric  transducer.  The  displacements  were  calibrated  using  a  Mi- 
chaelson  interferometer.  The  CT's  response  to  the  resonating  rod  was  then  measured  to  obtain  a 
calibration  of  the  CT. 

frequency  (to  about  1  MHz,  our  range  of  interest  for  this  experiment).  A  more 
complete  calibration  of  the  CT  has  been  completed  and  is  forthcoming. 

Data  from  the  CT  is  recorded  on  a  digital  oscilloscope  and  subsequently 
transferred  to  a  workstation  class  computer.  Time-signal  averaging  as  well  as 
discrete  Fourier  transform  analysis  is  used  in  the  interpretation  of  the  data. 


The  Samples 

All  samples  consist  of  standard  window  pane,  float  glass,  with  crossections  of 
25.4  mm  high  by  5.6  mm  wide  and  varying  lengths.  The  glass  has  the  following 
properties:  compressional  wave  velocity  of  5.64  mm/jus,  shear  wave  velocity  of 
3.35  mm/gs,  and  density  of  2610  kg/m3  which  corresponds  to  a  Young's  modu¬ 
lus  of  72  GPa,  and  a  Poisson’s  ratio  of  0.228. 

The  manner  in  which  the  largest  sample  (2.28  m  long)  is  supported  is 
illustrated  in  Fig.  1.  Urethane  foam  is  used  to  acoustically  isolate  the  sample 
from  the  aluminum  support  structure.  After  initial  experiments  to  test  the  ex¬ 
perimental  set-up,  a  vertical  cut  6  mm  deep  and  0.5  mm  wide  was  made  in  the 
sample  1.4  m  from  one  end,  effectively  creating  a  configuration  suitable  to 
detect  the  surface  breaking  crack  (see  Fig.  7).  Shorter  samples  were  used  and 
supported  in  a  manner  which  allowed  access  to  both  the  top  and  bottom  of  the 
sample. 


Experimental  Results 

In  support  of  the  theoretical  treatment.  The  two  dimensionality  assumption 
was  validated  through  moving  the  point  source  across  the  width  of  the  sample 
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Fig.  5.  Confirmation  of  two-dimensional  plate  assumption.  The  impact  point  of  the  steel  bearings  is 
moved  across  the  5.6  mm  sample  in  I  mm  increments.  No  appreciable  difference  was  observed 
unless  the  balls  were  glancing  off  the  edge  of  the  sample. 


Experimental 


Time  [ms] 


Theoretical 


5H 


Time  [ms] 


Fig.  6.  Experimental  data  and  theoretical  simulation  of  the  surface  displacement  at  5H  and  8H 
without  a  crack.  The  experimental  traces  are  multiplied  with  a  cosine  window  to  ease  theoretical 
computation.  The  experimental  surface  response  at  5H  was  used  with  the  Green's  function  to 
compute  the  source  function.  This  source  function  was  then  used  for  modeling  the  surface  displace¬ 
ment  at  8H.  The  analogous  technique  was  used  to  calculate  the  surface  displacement  at  5H. 


and  observing  the  signal  resulting  from  these  impacts.  Fig.  5  shows  the  experi¬ 
mental  arrangement  as  well  as  the  observed  signals.  After  propagating  through 
the  thickness  of  the  sample  the  signal  has  two-dimensional  character,  i.e..  it  is 
independent  of  where  it  originates  along  the  width  of  the  sample. 
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Fig.  7.  (a)  Experimental  set-up  to  study 
the  scattering  in  the  vicinity  of  the  crack. 
Receiver  locations  are  positioned  at  1  mm 
increments  (H  =  24.5  mm).  Experimental 
results  are  shown  in  Figs.  8-13.  (b)  Theo¬ 
retical  calculations  in  the  vicinity  of  the 
crack  were  made  using  the  finite  element 
mesh  in  region  R,  of  Fig.  3. 


(b) 


The  nature  of  the  source,  the  impact  of  a  steel  ball  bearing,  is  inherently 
simple.  Its  shape  under  the  two-dimensional  assumption  was  determined  from 
signals  recorded  at  5H  and  8H  from  the  source  (H  is  the  thickness  of  the  sample 
25.4  mm,  see  Fig.  7).  These  experimentally  observed  signals  together  with  their 
theoretical  counterparts  are  shown  in  Fig.  6.  The  experimental  results  at  5H 
were  used  with  the  Green’s  function  to  compute  the  source  function.  This 
source  function  was  then  used  for  modeling  the  surface  displacement  at  8H. 
The  same  technique  was  used  to  calculate  the  theoretical  surface  displacement 
at  5H. 

Scattering  in  the  far  field  was  examined  with  the  arrangement  shown  in  Fig. 
7.  The  source  was  located  at  a  distance  5H  from  the  crack.  The  receiver  was 
positioned  at  5H  -  5  mm  to  5H  +  5  mm  in  1  mm  intervals.  The  measured 
surface  displacements  for  the  cases  with  the  crack  and  without  the  crack  are 
shown  in  Fig.  8.  The  surface  displacement  changes  only  very  little  over  this 
range  if  no  crack  is  present.  The  forward  scattering  is  the  most  notable  effect  of 
the  crack  (receiver  locations  6-10).  Amplitude  spectra  were  calculated  for  each 
trace  to  allow  comparison  with  the  theoretical  calculations.  Figure  9  shows  the 
spectra  of  the  surface  displacements  at  the  positions  5H  +  1  mm  and  5H  -  1 
mm  with  and  without  the  crack  respectively.  The  experimental  data  as  well  as 
the  theoretical  results  are  shown.  The  experimental  spectra  were  normalized  to 
their  maximum  value  and  the  theoretical  values  were  scaled  in  such  a  way  that 
the  value  at  44  kHz  is  identical  with  the  experimental  value.  Fig.  10  shows  the 
same  spectra  in  the  range  100  kHz  to  200  kHz.  It  is  this  frequency  range  where 
the  introduction  of  the  crack  causes  the  greatest  changes  in  the  spectra. 
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rig.  s.  Normal  surface  displacements 
as  a  function  of  time  with  and  with¬ 
out  a  crack.  The  displacements  were 
measured  over  5  mm  on  either  side 
of  the  crack  at  1  mm  intervals.  The 
source  was  5H  (H  =  25.4  mm)  from 
the  crack  position.  The  upper  figure 
shows  the  results  without  a  crack  and 
the  lower  figure  with  a  crack  (See 
Fig.  7.  for  experimental  set-up). 


The  ratio  of  the  spectra  at  the  positions  5H  -  1  mm  and  5H  +  1  mm  with 
and  without  the  crack  are  shown  in  Fig.  1 1 .  In  Fig.  [2  the  ratio  of  the  spectra  for 
the  no-crack  and  the  crack  case  are  given  for  the  positions  5H  -  1  mm  and 
5H  +  1  mm.  Both  representations  show  a  distinctive  change  in  the  spectrum  in 
the  frequency  range  100  kHz  to  300  kHz.  In  Fig.  13  experimentally  obtained 
spectral  amplitudes  are  plotted  as  a  function  of  receiver  position.  Each  trace  is 
normalized  to  its  maximum  value.  The  presence  of  the  crack  introduces  a 
distinctive  maximum  of  the  spectral  amplitudes  in  certain  frequency  ranges. 

In  search  of  a  crack.  To  completement  the  two-dimensional  far  field  studies 
reported  in  the  previous  section,  we  performed  two  additional  experiments.  In 

- ► 

Fig.  9.  Amplitude  spectra  of  normal  surface  displacements  1  mm  on  either  side  of  5H  from  the 
source  with  and  without  a  crack.  5H  is  the  distance  between  the  source  and  the  crack.  Movement 
of  the  receiver  location  by  2  mm  changes  the  amplitude  spectrum  very  little  without  a  crack 
(bottom  figures)  and  an  appreciable  amount  with  a  crack  (top  figures).  The  solid  lines  are  experi¬ 
mental  data,  the  crosses  are  the  theoretical  results. 

Fig.  10.  Amplitude  spectra  of  Fig.  9.  shown  expanded  between  100  kHz  and  200  kHz  where  the 
crack  causes  the  greatest  changes  in  the  spectra. 
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5H  -I  mm  From  Source/Crack  5H+I  mm  From  Source/Crock 


5H- 1  mm  From  Source /Crack  5H+  Imm  From  Source  /  Crock 


5H-lmm  From  Source/NoCrack  5H  +  lmm  From  Source/NoCrock 


See  previous  page  for  figure  captions. 
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Fig.  11.  Ratio  of  the  spectra  at  position  5H  -  1  mm  and  position  5H  +  1  mm  without  a  crack  and 
with  the  crack. 
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Fig.  12.  Ratio  of  the  spectra  with  a  crack  and  without  a  crack  at  position  5H  -  1  mm  and  position 
5H  +  I  mm. 


one  case  the  source  was  stationary  and  the  receiver  was  moved  in  large  inter¬ 
vals  (40  mm,  nearly  2H)  across  the  sample.  In  the  data  shown  in  Fig.  14  the 
ultimate  crack,  i.e.,  the  end  of  the  sample,  is  clearly  visible.  No  effect  of  the 
induced  crack  is  discernable  however.  We  confirm  that  indeed  to  find  small 
cracks  one  must  look  in  detail  on  a  scale  that  is  on  the  order  of  the  crack 
whether  that  is  in  the  far  field  as  described  above  or  in  the  near  field  as  we 
describe  now. 

For  the  near  field  experiments  we  placed  the  receiver  and  the  source  on 
opposite  sides  of  the  sample  in  order  to  observe  the  near  field.  The  base  plate  of 
our  transducer  does  not  allow  the  source  to  be  closer  than  approximately  50 
mm  from  the  receiving  probe  when  both  receiver  and  source  are  on  the  same 
surface.  Figure  15  shows  the  recorded  surface  displacements  if  no  crack  is 
present  over  a  distance  of  180  mm  from  the  source  position. 

In  Fig.  16  we  show  the  experimental  set-up  which  we  used  to  search  for  the 
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Fig.  13.  Amplitude  variation  of  several  frequency  groups  as  a  function  of  position  on  the  sample. 
The  spectra  were  produced  from  the  traces  shown  in  Fig.  8. 


Crack  End  reflection 

TTti"i  r  T'T  T-r-rri  TTT-,-,  ^ -r  r  v  *  «  '  1  '  '  '  ▼I 


T=T=T=1  84  7  3 
26  33  41  49  57  66  74  82  1 

Time  [ms] 


Fig.  14.  Normal  surface  displace¬ 
ments  as  a  function  of  time  with  and 
without  a  crack.  The  normal  surface 
movements  were  measured  at  40  mm 
intervals  without  a  crack  (top)  and 
with  a  crack  (bottom).  The  distance 
between  the  first  receiver  position 
and  the  source  was  5H  in  both  cases. 
In  the  lower  figure  the  crack  was 
located  I  mm  beyond  the  first  mea¬ 
surement  location.  While  the  reflec¬ 
tion  from  the  end  of  the  sample  is 
recognizable  no  difference  in  the 
waveforms  due  to  the  crack  is  per¬ 
ceptible. 


crack  in  the  near  field.  Fig.  17  shows  the  two  sets  of  four  traces  recorded  for  the 
cases  with  and  without  a  crack.  The  differences  between  the  two  data  sets  are 
striking.  The  crack  is  easily  detectable.  In  traces  2  and  3  and  to  a  lesser  extent 
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Fig.  IS.  Near  field  experi¬ 
ment  in  the  no-crack  case. 
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Fig.  16.  Experimental  configuration  used  to  study  the  effect  of  the  crack.  Normal  surface  displace¬ 
ments  on  the  bottom  of  the  sample  are  measured  in  response  to  ball  bearing  impacts  at  locations 
1,2,3,  and  4  above  the  sample.  The  results  are  illustrated  in  Fig.  17. 


in  trace  4  the  presence  of  the  crack  prevents  the  destructive  interference  which 
occurs  during  the  first  0.03  ms  in  its  absence.  For  the  no-crack  case  the  data  in 
Fig.  17  represents  a  subset  of  those  shown  in  Fig.  15.  There  too.  the  regions  are 
clearly  recognizable  over  which  the  destructive  interferences  occur,  i.e.,  at  the 
20  mm  and  25  mm  traces  and  again  at  the  45  mm  and  50  mm  traces. 
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Fig.  17.  Normal  surface  displacements  on  the  opposing  side  of  the  sample  from  the  source.  The 
source  is  placed  at  one  of  the  four  positions  shown  in  Fig.  16.  The  receiver  is  stationary  and  on  the 
opposite  side  of  the  sample.  Times  after  0. 12  ms  are  subject  to  reflections  from  the  ends,  and  should 
be  ignored. 


Conclusions 

In  this  study  we  have  presented  experimental  and  theoretical  treatment  of  a 
two-dimensional  wave  propagation  problem  in  a  plate  with  and  without  a  sur¬ 
face  breaking  crack.  For  the  crack  free  case  we  showed  good  correspondence 
between  the  experimentally  observed  normal  displacement  time  series  and  the 
corresponding  theoretically  calculated  ones  when  the  reciever  was  located  in 
the  far  field  of  the  source.  Finite  element  based  calculations  of  changes  in 
amplitude  spectra  in  the  vicinity  of  a  crack  are  generally  in  good  agreement 
with  experimental  observations. 

High  detectability  of  the  crack  was  achieved  when  the  source  and  receiver 
were  located  in  the  near  field  of  the  crack.  While  these  observations  were  not 
treated  theoretically  in  this  study,  together  with  the  theoretical  work  they  may 
offer  the  basis  for  the  design  of  a  practical  crack  detection  gauge.  Both  the 
source  and  the  receiver  could  be  housed  in  the  same  physical  unit  and  thus  be 
easily  transportable  and  usable. 
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Appendix  I 

Let  us  write 


a  -  k\ m  -  2k2,  b  =  2kr)tm,  c  =  2kr)2m 
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Then  the  matrix  (/mU)]  have  the  elements 

IT,  =  ike111'"1.  In  =  -ivin,eiV2"z 
IT i  =  ike,1""{h  l\  ITa  =  it]2me,r)1"{k~l) 

IT,  =  In  =  ikeir)lml 

lii  =  -iVime'7"Jh~z\  ITa  =  ikeir,2"(h~z> 

IT,  =  -a/xmeir""z,  lyi  =  -C/J.meir)1"z 

IT  =  -aixmeill'Jh-zK  ITa  = 

IT,  =  b /J-meITlt"z ,  1 42  =  ap.me,r>-"z 
IT  =  bnmeil>'"'h~z),  ITa  =  atimeil<'"^z) 

The  elements  of  the  4  x  1  matrix  [  Tf]  depend  on  the  direction  and  location  of 
the  source.  If  the  source  is  in  the  jc-direction  at  (*',  z')  in  the  mlh  layer  then, 

YT\  =  ikA,meiJI'Jz'z 1  -  ir)2mB>meiri-"iz~z'1 
YT,  =  ±ir)lmAin,eiTI'"lz~z'}  ±  ikB[me‘T'-Jz'z'1 

—  YT  =  -aA^e*1""1^"1  - 

1  _  yjl  ~ 

On  the  other  hand,  for  the  force  in  the  ^-direction  we  have, 

YT  =  ±ikAime'1,,mk~z^  ■+  irnmB^me,r,'-"z~ z 
YT  =  i'n,mAimeil,'"k~zl  +  ikBime,r,'"]z~z' 

_L  yt  =  +aAimein>^:~z^  —  cBimeil,'mk~z 1 

l  YT  =  _ 

Hm  4 
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